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^ ' Abstract 

00 ! 

\^ \ We give a complete analytical computation of three-point one-loop integrals with one heavy 

' propagator, up to the third tensor rank, for arbitrary values of external momenta and masses. 

(N : 
o . 

■ 1 Introduction 
O ; 

' The study of the dynamics and spectroscopy of hadrons containing a heavy quark has been greatly 

simplified and systematized with the introduction of heavy quark effective theory (HQET) [p. 
Qh. Heavy-particle theories along similar lines have also been successfully applied in other, related 

(— I I contexts. Thus, in those cases where a chiral approach to the strong interactions of heavy hadrons 

■ with light mesons is applicable, a combination of chiral and heavy-quark symmetries leads to 
•rH , heavy hadron chiral perturbation theory (HHChPT) A heavy-particle expansion has also been 

I developed in the chiral-perturbative framework for nucleon-meson interactions, which constitutes 

' the so-called heavy baryon chiral perturbation theory (HBChPT) |2|. 

In this paper we report on a complete analytic computation of three-point loop integrals with 
one heavy propagator, up to the third tensor rank, for arbitrary real values of masses and residual 
momenta. Together with the results for vertex integrals with two heavy propagators |^] , this gives 
the set of all one-loop three-point integrals in heavy particle theories. 

For the calculations we use the same method as in , namely, we compute the vertex integrals 
with heavy propagators as limits of ordinary three-point integrals. The latter are evaluated by 
closely following the method of Q, so that the calculation is standardized and the results expressed 
in terms of the usual logs and dilogs. Since we only have to compute to leading order in a certain 
limit, however, the calculations are drastically simplified as compared to ordinary vertex integrals. 

In the next section, we define the integrals to be studied and give the calculational details and 
results for the scalar integral. We also discuss some particular cases which are of importance in 
applications and as cross-checks of our results. Then, in section |3| we consider the vector, second-, 
and third-rank tensor integrals, which are given in terms of scalar integrals. 

* E-mail: abouzas@mda.cinvestav.mx 
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2 Vertex integrals 



The three-point loop integrals we consider are of the form, 

^ (27r)'^y {2vi-5M + ie)i{i-pf-m^ + ie){{i-p')^-m'^ + ie)' 

Here is the velocity of the heavy lines, p'^ and p''^ the residual momenta of the external heavy 
particles, and 6M the residual mass of the internal one, given by the mass splitting within the 
corresponding heavy-quark symmetry multiplet. m and m' are the masses of the light particles 
in the loop, which are pseudoscalar mesons in HHChPT and HBChPT, and (massless) gluons in 
HQET. The scalar integral is finite, whereas the vector and higher-rank tensors are divergent, with 
degree of divergence n + d — 5. We consider tensor integrals with n = 1,2,3 below. These are 
computed in terms of integrals with lower rank and fewer points. 

By shifting the integration variable i ^ £ + {p' + p)/2 , the scalar integral can be written as. 



'^^ {2t:Y j'^ \2vl-A + i£) {{I + qY-m? + ie) {{I- qY-m'-^ + iey 

with A = 5M — V ■ {p' +p), q = {p' —p)/2. Thus, Ti.3 depends explicitly on g^, through q^ and v ■ q, 
but depends on (p' +p)^ only through A. Since TC^ is invariant under the transformation q —q, 
m <-> m' of its arguments, we can choose v ■ q > without loss of generality. 

We will compute the scalar integral TC3 as the large-mass limit of the auxiliary integral 

" (2^ I'^^Uie + Mt>)2 - (M + A/2)2 + ie) {(i + qf - m? + is) {{£ - qf - m'2 + ie) ' 

Since both TL^ and ^3 are convergent, we clearly have, TL^ = \\m.M-*oo J^I'^-^. The calculation of 
to leading order in M will be carried out by closely following the standard method of Q. We 

will only sketch the main steps of the calculation. 

Starting from the usual Feynman parametrization for Ti^, after integrating over and over 

the Feynman parameter multiplying the last propagator in (|^, we obtain. 
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D = M'^z^ + Mz(-2(1 - 2x)v ■ q + A) + xm^ + (1 - x)m'^ - 4x(l - - ie. 

In this expression we have retained only the leading terms in M in the coefficients of the polynomial 
in the denominator. Clearly, only the region where z = 0{M-^) wiU contribute to Tis to 0{M-^), 
so we have also set the upper limit of the integral over x to 1 instead of 1 — z. 

We notice that if we make the change of variable A = Mz in (^) and let M ^ 00, we are led to 
the HQET parametrization for Ti^, confirming the validity of the approximations leading to @. 

Following Q, we make the change of variable z ^ z — ax/M, with a one of the roots to 
{av - 2qf = 0, 

a±=2{vq±H). (5) 



We henceforth assume that cf' < {y ■ q)^, so that |q| = V {v ■ q)"^ — cp' > is real. For the change of 
variable we choose a = a+ > 0, leading to, 

~ (47r)2 io '^^ J ax/M M2z2 _ (4M|q|z - zo)x + Mz{-2v • g + A) + m'2 - ' 



2 



with zq 



-{m 



— m^) — a^_(A — 2|q|). Exchanging the order of integration, 



1 



(47r)2 

X 



a/M 



dz 



Mz/a 



dx + 



a/M 



dz I dx + 



1+q/M 



dz 



dx 



(7) 



^2^2 _ (4M|q|2; - zo)x + Mz{-2v ■ q + A) + m 



/2 



le 



In the last integral in this equation z = 0{1) over the entire domain of integration, so that the 
integral itself is C(M~2) and can be neglected. We will drop it henceforth. Integrating over dx we 
then have, 



X 
Y 
Z 



1 



(4vr) 



a/M 



dz- 



-1 



4M|q|2: - Zq 



log ( ^ ) + 



1 



dz- 



-1 



M^z^ _^ Mz{-2v ■q + A) + m' 



12 



12 



M^z^ + Mz{-2v ■q + A)+m' 
M'^z^ + Mz{-2v ■q + A)+m'^ 



ie 
ie 



(47r)2 J^iM 4M|q|2: - zq 
(4M|q|z - zo) 
M 



log(§ 



(8) 



a 



-z(4M|q|z - Zq) 



We notice that there is no singularity at the zero of the denominator of the integrand in (|^, 
z = zo/(4M|q|), since X = Y = Z there and the logs vanish. The roots of these polynomials in z 
are given by. 



Xl,2 



yi,2 



Zl,2 



V ■ p' + 2\q\ 



6M 



V ■ p ■ 



+ ie 



— ^4g^ + m'^ — lb ^ ^4q2 — (771' + m)^^ ^4g2 — (^rri' — m)2 j + iea^ 



2a 



V ■ p 



(9) 



5M 



V ■ p' 



JM\2 



m'2 + ie, 



where in the expression for yj we denoted a = sgn{q'^). The fact that X, Y, and Z are all equal at 
z = Zo/(4M|q|) leads to the important identity. 



4|q| 



Xl 



4|q| 



X2 



Zq 
4|q| 



yi 



Zq 
4|q| 



y2 



Zq 
4|q| 



zi 



Zq 
4|q| 



Z2 



(10) 



Factorizing the polynomials in the arguments of the logarithms in (H), we get 



1 



+ 



(47r)2 
1 



a/M 



dz 



dz 



-1 



4M|q|z - Zq 
-1 



(47r)2 J^/M 4M|q|z - zq 



log 
log 







yi 




M 


(- 







m 

M 



log 



Zl_ 

M 



X2\ 

Mj 



log 



MJ 



z 



z - 

Z2_ 

M 



Z2_ 

M 



■ (11) 



In order to be able to distribute the integrals inside the braces without introducing spurious 
singularities at the zero of the denominator, we add and subtract the value of each logarithm at 
z = zo/(4M|q|), with the help of identity (|10|). We can then use the fact that 2:1^2 have imaginary 
parts of opposite sign, and analogously yi^2 and 2:1^25 to split the logarithms. (In the case of the 
logarithm containing yi^2, proper account must be taken of the fact that there is a factor a^/a+ in 
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the argument. Whereas a+ > 0, we have a_ > if > 0, and a_ < if < 0.) This procedure, 
which is described in detail in |Q| (see also Q, and appendix^), leads to a set of integrals that can 
be evaluated in terms of dilogarithms, with the result, 
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^3 = Ji?oo^^3 = (i^o^Sd {J'i{y,)+J'2{x,)-T^{z,)) (12) 



i=i,2 



^3(.) - ^i(x) +^,(x) = -LiJ^^) - llog^ -^1^1^ 



zo — 4|q|x/ 2 V 

In this equation the functions J-2 and J-^ are defined in terms of an arbitrary mass parameter 
/i > 0, analogous to the dimensional regularization mass unit. Tij, does not depend on /i since, 
by virtue of the identity (|lO|), all // dependence cancels in the sum in ([T^). In principle, /x can be 
eliminated from 7^3, but we prefer not to do so. In appendix ^ we discuss the form of J-2 in more 
detail, and give several equivalent forms for it, some of which will be used below. (|l2|) is then our 
final result for 7i^. 

The cuts of Tiz as a function of external momenta can be easily found from ([Tl|). As M ^ oo, 
the integral of the logarithm containing yi^2 has a cut at (m' + m)^ < Aq^ < +oo, that of the log 
containing xi^2 has a cut at 5M/2 + m < v ■ p < +oo, and the integral of the log containing zi^2 
has a cut at 6M/2 -\-m' < v ■ p' < +oo. 

2.1 The massless case 

In HQET the heavy quark-gluon vertex at one loop is expressed in terms of integrals of the 
form dl]) with m' = = m, and 6M = if the heavy quark residual mass term in the Lagrangian 
is chosen to vanish. The scalar integral Ti^ has been given in Q in that case, for < 0, v ■ p' < 0, 
V ■ p < 0. It is then of interest to evaluate our result (|l^) in the massless case, both because of its 
applicability to HQET and to cross-check our calculation with that of [||. 

Defining w' = v • p' , w = v • p, Q = 2|q| = |p' — p| > 0, a^i = sgn{w') and aw = sgn(?x;), and 
setting m' = m = 6M = 0, from (^) we get, 

^3 = 77^7^ jLiaf ) - ? (^{w' + w + Q){w' - w - Q) - iea, 

(47r)^ 4|q| L \w' + w + Q + ieaJ 2 V/" 

- Li2 ( ] - I log' (^{w' -w + Q){w' + w-Q) + teaw] (13) 

\w' + w — Q — ieaJ 2 Va* / 

- Li2 1—, , , —TT. : + o log [ ^{w + w - Q)[w - w - Q) - leaw' 

\(w' - w - Q){w' + w + Q) - leaw J 2 V/" 

^ / {w' + W + Q){w' -W + Q) \ 1, 2/^1// I ^^ ■ 

+ Li2 -n\rr' 7^ ]+ -\og^ (—{w' + w + Q){w' -w + Q)+ieaw 

\{w' + w — Q){w' — w — Q) — ieaw> J 2 VA* 

This expression is in full numerical agreement with eq. (24) of ||5| in the region where g^, w and w' 
are all negative. Also, ([l^) has cuts at if' > 0, w > and > 0, as it should. 
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2.2 The limit g ^ 

In the limit q ^ the integral Ti^ can be easily computed as the difference of two-point integrals. 
Using the notation of Q (see appendix we have, 

H3\q=o = —p^ ^ (j2(A,m'^) - 2'2(A,m^ 

— V 

.2 

2 



1 



1 



(47r)2 m'2 — 



.7^(2;) = — 1 + ie \o<g{x — \f . 



- 1 + ie^ — log(^x + — 1 + 



(14) 



where ^' = A/(2m'), C = A/(2m). 

In order to obtain the value at g = of T^s directly from (^) we have to expand the functions 
in the sum to first order in ^ due to the factor l/|q| in front of the sum. 

We consider first the terms Z]j(-^i(?/j)~-^i('2^i)) ™ (O)' ^'^^ limit g — > the two dilogarithms 
in the definition of J-\ (see (p!2|)) cancel each other up to terms of 0(|qp), so these terms give a 
vanishing contribution to Hj, at g = 0. It should be noted, however, that the expansion about 
(7 = of yi^2 is singular because of the factor l/a_ in their definition (P). Thus, we expand first 
about |q| = 0, (7'^ > and take the limit g*^ — > afterwards. There are no such complications in 
the case of T\{^z\^2). 

In order to expand the remaining terms X]j(-^2(a^j) — •^2(-2j)) in (12) about g = 0, we use the 
form ( [A. 5 ) for T2{x). To lowest order, we get. 



4|q| 



^(.F2(x,)-^2(^,)) + 0(g) 



i=i,2 



i=i,2 





+ ^log 


(-1) 


I ^0 







(15) 



Here, the zeroth-order terms in the expansions of Xj, Zj, zq about ^ = 0, respectively. 

Using the notation of (|lj), we have. 



2:1,2 



m 



Zl,2 



m 



Zo 



+ (16) 



Substituting (16) into (15) we recover ([l^). 



3 Tensor integrals 

In this section we consider tensor three-point integrals, up to the third rank. These are given in 
terms of integrals with smaller ranks and fewer points, by using the well-known method of Q. The 
general tensor integral has the form (|l|) . We will restrict ourselves to integrals of the standard form 

j.A-d f pa\ . . . ooin 

3 ^""^ ^' ^' ^ (27r)rf J {2ve -A + ie) {{£ + qf - rrfi + ie) ((£ - qf - m'2 + ie) ' 

(17) 

By shifting t+{p' +p)/2 in (|l|) we can express Wg^""" in terms of Hg^""" as, 

n 

^ H;^^-""(z;,g;A,m,m') + 5]r^"i •••r"^H3^+' ■"">(7;,g;A,m,m'), (18) 

i=i 
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where = l/2{p' +pY, and y4{"i"2-a4 = ^aiaa-a. _^ ^a2-asai _^ _^ _ Clearly, the 

scalar integral Ti^ = H3. The standard form (|T7|) has the advantage that it depends explicitly on 
only two four- vectors, v^^ and g'^, instead of three as in (|l]). 

3.1 Vector integral 

For the vector integral we write H3 (f , q; A, mi, 7712) = Viu" + V2g'"/|q|, with |q| = V ("V • qY ~ 
and Vi^2 = Vi^2('y • 9, 9^, A, '"T'I, 'm.2)- If |q| = 0, then (7" oc and we can set V2 = 0. If |q| ^ 0, 
then, 

H^yi = -q\^H1 + vqq^H'^, |q|V2 = f -g 7;„H? - g„H?, (19) 

with, 

?;q,H3 A,mi,m2) = ^So(4g^, mi, 7712) + y H3(t>, g; A, mi, m2) 

gQ,H3(u,g; A,mi,m2) = ^J2(A - 2t; • q', m2) - ^X2(A + 2t! • g, mi) (20) 

m? — mo , , , . , 
H ^ H3(v, g; A, mi, m2). 

The scalar two-point integrals X2 and i?o ^-re given in the appendix. 



3.2 Second-rank tensor integral 

The decomposition of Hg'^ in terms of form factors can be written as, 

Hf{v, q; A, mi , m2) = Tig"'^ + j^y^^^P + ^3^ + T^—^. (21) 

The components of Hg^ are obtained by direct computation, 

Fi = 5a/3Hf = ^UA -2vq, ma) + ^T2(A + 2vq, mi) + "^1 + ^1 - V 

F2 = VaVfsH'^'^ = -V ■ qBi (^Aq'^,mi,m2^ + ^ (^y ~ ^"^^ -^0 (4g^,mi,m2) + (22) 

F3 = ^a^/jHg^ = Y (2^2'('"> A - 2w • g', m2) - 2:^('V, A + 21" • q', mi) 

2 22 
+ ^ (T2(A -2vq, m2) + T2(A + 2^; • g, mi)) + ""^""'^ g^Hg 

^2 ^ 

F4 = ^^ag/3H3^ = -g^5i(V,mi,m2) - ySo(V,mi,m2) + —qa^s 

where, on both sides of these equations the arguments (f , g; A, mi, m2) of H3, H3 and Hg^ have 
been omitted for brevity. To one-loop level, the form factors in (^Tj) are given in terms of the Fs 

as. 



2|q|2Ti = (1 + (|q|'Fi + q^f2 + F3 - 2vqfi 



\4 1/ / \\ 2 

e\ w ^ 1 /„/ \2 ' - - sQ v-( 



2|q|^T2 = I + -U-F1+ 3 + -^F2 + ^ 2(.;.g)^+ 1 + ^U^ F3-(6 + 6)VifF, 



27 Iqr" |q|2 V ' V 2; V |ql2 



4 



6 



2|qpT3 
2|q|'T4 



1 + o lql Fi + 2{v.qY + 1 + F2 + 3 + 



(6 + e)v-q¥4^ 



(23) 



l + -h;-g|q|Fi- 3 + 



3+^)fTF3 + A(V + (4 + 6)(..,)2) F4. 
2/ |q| |q| ^ ^ 



Equations (^) and (23) give an explicit expression for H 



al3 



3.3 Third-rank tensor integral 

The decomposition in terms of form factors can be written as, 
H3^''(w,g; A,mi,m2) 

These form-factors are given impUcitly by the following linear relations, 

{d + 2)Si + {d + 2)vqS2 + SvqSs + {q^ + 2{vqf)S4 + S5 + q^v-qSe = gx,,v^y^3 
{d + 2)vqSi + {d + 2)q^S2 + {q^ + 2(vqf)S-i + Sq^v-qS^ + vgSg + q^Se 
SvqSi + (g2 + 2{vqf)S2 + (q^ + 2{vqf)S3 + (2g2 + {vqf)vqS4 

+v-qS5 + q'^v-qSe = vxv^qu^s 



(24) 



Xfiu 



Xiiu 



(25) 



(g" + 2ivqy)Si + 3g"(^;-g)S2 + v-q{{vqy + 2(7^)83 + <7'(r + 2{v-qy)Si 



Xfiw 



+{v- qyS5 + q'^v-qSe = vxqf^Qu^ 

3Si + 3v-qS2 + 3^;-g'S3 + 3(t;-q')^S4 + S5 + {v-q^Se = vxVf_,Vi,^^'^'^ 
•iq^vqSi + 3/S2 + ?,{v-qfq^S^ + 2,[vq)q^Si + [v-qfS^ + g^Se 

where the right-hand side can be obtained by direct computation. 



gx^Qu^s^" = ^{inv,A + 2vq,mi)+I^{v,A-2vq,m2)) 

+ |- (J2(A - 2t;-(7, m2) - + 2vq, mi)) + "^1 + ^2 " 2g 



1 



A 



2 B^'^{2q,mi,m2) +2q vqBi + -v-qq -Bq + yi^a^/jH 



0/3 
3 



^q^qpE'^'^ {2q,mi,m2) + 2g'^Si + ygag/jHg^ 



2 



(26) 



-VaVf3B2'^{2q,mi,m2) + -{v-q) Bo + 2{vq) Bi + —VaVf3H^ 



jQaqf3 {^2^{v, A - 2vq, 7712) - 12'^ {v, A 2wg, 



mi 



+^g„(2'2"(^^,A-2t;-g,m2) + T2"(^^, A + 2t;-g, mi)) 

4 22 
+ ^ (J2(A - 2vq, m2) - 2'2(A + 2vq, mi)) + ~ ""^ gag/3 H 



a/3 
3 ■ 



In these equations i?o = BQ^Aq'^ ,mi,m2) and similarly i?i, and H3 = H3(?;, g; A, mi, m2) and 
similarly H? and H?''. 
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4 Final remarks 



In phenomenological applications, the exact functional dependence of Feynman integrals on masses 
and residual momenta is usually not needed. Often, the first few terms in a series expansion in 
some of the parameters provides the required accuracy. We believe, however, that the exact 
analytic computation presented here does not require more calculational effort than approximate 
schemes. It has the added advantage of being valid over the entire physical region for internal 
and external masses. Applications of the results presented here to loop graphs in heavy baryon 
chiral perturbation theory which have a calculable dependence on the ratio of the A-nucleon mass 
difference to the pion mass will be given elsewhere. 

Our approach, which is based on well-known methods and results for vertex integrals in renor- 
malizable theories [^, ^, streamlines the computation so that it can be easily reproduced and 
verified, and adapted to other, more complicated one-loop diagrams. Another important feature 
of this approach is that contact with the unitary cuts of the diagram, as given by Cutkosky rules, 
is explicitly kept at every step of the calculation Q . 
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A Calculation of 



In this appendix we discuss the calculation of the second integral in (11). We consider first the 
basic integral, whose computation is standard 

dz 



a/M 4M|q|2; 




where J^2{x) is defined by this equation. In (|A.l|) we have made use of the relation Li2(x) = 
— Li2(l/x) — 7r^/6 — l/21og^(— x) in order to obtain the leading M dependence of Li2(4M|q|/(zo — 
4|q|x)). 

Denoting by Q the second integral in (0), we have, 

(A.2) 



Clearly, those terms in J-2{x) which do not depend on x will cancel in ( [A.2| ), so they can be omitted. 
Also, due to the identity ([lO|), Q does not depend on the value of M, which can be replaced by an 
arbitrary mass /x > that remains constant as M ^ oo. Thus, in ( A. 21) we can write. 



zo — 4|q|x / Vzo — 4|q|x/ \ ji J V-^o — 4|q|x/ V 4|q|/x 

2 U-4|q|i/ '^\z„-i\<i\x) '^\7.a-iHx) ^ ' 

The last two logarithms in this equation can be rewritten according to the identity log^(— z) — 
21og(— z) log(z) = — TT^ — log'^(2;), which is valid for the principal determination of the log. Dropping 
the constant term, we see that in ( [A.2 ) we can write, 

(x) = Li2 ( + log ( log f ^) - log ( log ^ ^0 - 41^1- 



zo — 4|q|x / Vzo — 4|q|x/ V /U / V-^o — 4|q|x/ V 4|q|/i 

hoe(^^). (A.4) 



2 \zo — 4|q|a;, 

This equation can still be rewritten by completing the last two terms to the square of a sum of 
logs. Once again discarding x-independent terms that do not contribute to (|A.2|), we get, 



^2(x) = Li2 ( M-ZA) + logf^^^^^) logf ^) + i log^ ( ^^^^f^ ) . (A.5) 



zq — 4|q|a; / V-^o — 4|q|x/ \ ^ J 2 Vzo — 4|q|a 

From (|A.4|) , by means of the identity 



zq — 4|q|a; / \ zq — 4|q|2; 



- - Li2 C-^^^] - logf^""'!^!"^ logf^^^^ 



Li2 -7—^ - log jT-r- log , , I I 

D V^o — 4|q|x/ V-20— 4|q|x/ V -^o — 4|q|x 



we obtain .7-2 (x) as given by (11), up to constant terms. All of these expressions for T2{x) lead to 
the same results when substituted into ([A.3|) or (111]). 
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B Loop integrals 

In this appendix we give a list of loop integrals used in the foregoing. More complete calculations 
can be found, e.g., in |^, ^, |^, ^ and references therein. Divergent integrals are separated in a 
dimensional-regularization pole term and a finite remainder. Ji = ij,V^7re~^. 



{2ttY J i2_^2^i^ 87r2e Idn'^ 
ao(m^) = log ^ - 1 



{2ttY J {l'^ -ml + ie){{l + pf -ml + ie) 



Svr^e 167r2 



[ , , ( , .ni\ ml , ,p \ 
/ axlog (l-x)^;5- + j;^;5 — x{l — x)^—ie\ 
Jo \ A*^ At^ /""^ / 



5f(p,mi,m2) = -— ^ / d I— 2 , ■ \ i(B , \2 2— T =P^Bi(p , mi, 7712) 

[zttY J \" — mf + le) {[t + pY — m2 + le) 

p^Bi{p^,mi,m2) = ^ (Ao(mi) - ^o("i2) - (p^ + - m|)So(p^,mi, 777.2)) 



B^^{p,m„m2) = j^,Jd ^ (^2_^2^,,) i^i + pY-ml + 7e) = ^ " 7^ j + "^1^ 

{d-l)G = ^Ao{m2) + mlBo{p'^,mi,m2) + ^{mj-ml+p'^)Bi{p'^,mi,m2) 

1 1 
L = -Ao{m2) - -{mj-ml+p'^)Bi{p'^,mi,m2) 



{27rY J i2ve- A + i£){i^ -m^ + ie) 



I^(v,A,m) = / d^i- 5 - = F(A,m)v 

F{A,m) = ^Ao{m) + ^l2{A,m) 



Tr(7;,A,m) = / d'^l . ^ ■ w,2 = loiAMg'"" + h{A,m)v>^v^ 

A . _ / o A2^ 



(d- l)/o(A,m) = --Ao(m)+ (^m^-— j J2(A,m 
(d- l)/i(A,m) = ^Ao(m) - ( - ^ ) J2(A,m) 
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